home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / BESSK.DEM < prev    next >
Text File  |  1991-04-29  |  652b  |  33 lines
  1. PROGRAM d6r25(input,output,dfile);
  2. (* driver for routine BESSK *)
  3. VAR
  4.    n,i,nval : integer;
  5.    val,x : real;
  6.    txt : string[27];
  7.    dfile : text;
  8.  
  9. (*$I MODFILE.PAS *)
  10. (*$I BESSI0.PAS *)
  11.  
  12. (*$I BESSI1.PAS *)
  13.  
  14. (*$I BESSK0.PAS *)
  15.  
  16. (*$I BESSK1.PAS *)
  17.  
  18. (*$I BESSK.PAS *)
  19.  
  20. BEGIN
  21.    glopen(dfile,'fncval.dat');
  22.    REPEAT readln(dfile,txt) UNTIL (txt = 'Modified Bessel Function Kn');
  23.    readln(dfile,nval);
  24.    writeln(txt);
  25.    writeln;
  26.    writeln('x':5,'actual':14,'bessk(n,x)':18);
  27.    FOR i := 1 to nval DO BEGIN
  28.       readln(dfile,n,x,val);
  29.       writeln(x:6:2,'   ',val:13,'   ',bessk(n,x):13)
  30.    END;
  31.    close(dfile)
  32. END.
  33.